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Abstract 

A  large  block  of  material  is  subjected  to  a  constant  average  strain  rate  resulting  in  deformation  in 
the  antiplane  shear  mode.  The  material  is  assumed  to  be  isotropie,  incompressible  and  hyperelastic.  To 
simulate  the  situation  in  which  strain  hardening  of  die  material  is  overwhelmed  by  thermal  softening,  the 
shear  stress  as  a  function  of  shear  strain  has  a  local  maximum.  The  block  contains  a  material  imperfection 
in  the  form  of  a  small  region  in  which  the  material  is  less  stiff  than  the  otherwise  homogeneous  block,  and 
in  which  the  maximum  possible  stress  is  less  than  in  the  surrounding  material.  The  transient  defmmation 
is  analysed  numerically  by  means  of  the  Unite  element  method.  It  is  found  that  a  shear  band  initiates  in 
the  region  of  the  material  defect,  and  then  propagates  in  a  crack>like  fashioa  into  the  rest  of  the  body.  A 
prediction  of  the  speed  of  propagation  of  the  edge  of  the  band  is  made  on  the  basis  of  a  linearization  of 
the  governing  field  equations  around  a  state  corresponding  to  the  strain  level  at  which  the  local  maximum 
in  stress  occurs. 
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1.  INTRODUCTION 

The  localizatkm  of  deformation  into  shem  bands  is  commonly  observed  during  high  strain  rate 
inelastic  deformation  of  metals  and  other  materials.  Such  bands  persist,  once  they  are  formed,  and  the 
subsequent  deformation  proceeds  in  a  markedly  nonnniform  manner.  Aside  from  their  intrinsic  scientific 
interest,  such  shear  bands  are  of  technical  interest  because  they  provide  a  deformation  mode  with  an 
unacceptably  low  level  of  energy  absorption  for  some  applications  and  they  are  often  precursory  to 
fracture. 

Experimental  studies  on  shear  bands  have  been  carried  out  by  Nadai  (1031),  Zener  and  HoUomon 
(1944),  Rogers  (1970)  and  more  recently  by  Costin  et  al.  (1070)  and  Moss  (1080).  Nadai  has  assembled  a 
series  of  photographs  of  Luder’s  yielding  in  steel  specimens,  which  demonstrated  that  the  intense  shear 
entails  the  crack-like  propagation  of  bands  in  highly  nonuniform  conditions.  The  high  strain  rate  and 
high  temperature  within  the  band,  as  reported  by  Rogers,  Costin  et  al.  and  Moss,  indicated  that  the 
deformation  is  nearly  adiabatic  and  that  a  phase  transformation  may  occur  during  deformation.  The 
high  temperature  is  thought  to  be  due  to  the  heat  generated  through  dissipation  of  the  mechanical 
work  done  to  deform  a  material  inelastically.  If  the  deformation  is  very  rapid,  the  heating  process  will 
be  nearly  adiabatic.  Local  heating  tends  to  reduce  the  resistance  of  a  material  to  further  deformation 
which,  in  torn,  leads  to  the  generation  of  more  heat,  and  so  on.  The  material  becomes  unstable,  in  the 
sense  of  mechanical  equilibrium,  and  a  highly  localised  intense  deformation  develops,  typically  in  the 
form  of  a  shear  band  aligned  with  the  direction  of  maximum  shear  stress. 

Analytical  studies  were  initially  reported  by  Thomas  (1961),  Hill  (1062)  and  Mandel  (1966).  In  the 
latter  two  cases,  shear  bands  were  shown  to  correspond  to  statbnary  acceleration  waves  in  elastic-plastic 
solids.  More  recently,  Rudnkki  and  Rke  (1075)  and  Rke  (1076)  addressed  this  problem  for  inelastk 
materials  by  seeking  conditions  for  whkh  the  constitutive  relations  could  allow  a  bifurcation  from  a 
homogeneous  deformation  state  into  a  concentrated  shear  band  mode.  In  a  series  of  articles,  Knowles  and 
Sternberg  (1978,  and  references  therein)  determined  that  the  formation  of  surfaces  with  discontinuous 
displacement  gradients  in  otherwise  homogeneous  nonlinear  elastk  solids  is  due  to  the  loss  of  ellipticity 
in  the  equations  of  continuing  equilibrium.  These  studies  were  concerned  with  the  possible  existence  of 
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localizations,  and  the  kinetics  of  formation  was  not  addressed.  An  analysis  of  the  evolution  of  a  shear 
band  under  rising  remote  loads  was  presented  by  Abeyaratne  and  TViantafyllidis  (1980).  They  introduced 
an  imperfection  in  a  large  block  of  nonlinear  elastic  material  (either  heperelastic  or  hypoelastic),  and  they 
examined  the  differences  between  the  uniform  deformation  field  under  rising  load  and  the  nonuniform 
field  due  to  the  imperfection  by  means  of  a  regular  perturbation  analysis.  They  concluded  that  the 
shear  localization  behavior  may  be  exhibited  if  the  incremental  equilibrium  equations  loose  ellipticity. 
Moreover,  they  found  that  the  effect  of  the  imperfection  on  the  deformation  field  was  very  small  until  the 
remote  uniform  stretch  reached  values  which  were  very  close  to  the  critical  value  at  which  the  equations 
for  uniform  tensile  deformation  loose  ellipticity. 

An  alternative  approach  to  the  study  of  the  criteria  for  the  onset  of  shear  localization  in  one 
dimensional  modeb  has  been  considered  by  some  authors,  including  Recht  (1964),  Argon  (1973),  Culver 
(1973)  and  Staker  (1981).  They  have  proposed  instability  criteria  (i.e.,  criteria  for  the  formation  of  a 
shear  band)  based  on  the  attainment  of  a  local  maximum  in  the  dependence  of  shear  stress  on  shear 
strain  during  a  homogeneous  deformation.  In  a  departure  from  the  notion  of  a  criterion  based  on  a 
stress  maximum,  Clifton  (1980),  Bai  (1980),  and  Burns  and  Trucano  (1982)  examined  the  growth  of 
infinitesimal  periodic  nonuniformities  in  an  oth«wise  uniform  simple  shearing  deformation  in  a  material 
exhibiting  strain  hardening,  thermal  softening  and  strain  rate  sensitivity.  Clifton  and  Bai  assumed  that 
the  basic  material  parameters  occurring  in  the  problem,  such  as  the  strain  hardening  and  thermal 
softening  parameters,  are  time  independent.  Thty  found  that  the  magnitude  of  the  nonuniformities  may 
grow  or  decay  in  time,  depending  on  the  material  parameters,  the  average  strain  rate,  and  the  fixed 
spatial  wavelength  of  the  initial  periodic  nonuniformity.  Burns  and  IVucano  extended  the  study  and 
allowed  the  material  properties  to  be  time  dependent.  They  found  that  the  stability  prediction  for  time- 
dependent  material  properties  is  somewhat  different  from  that  based  on  time  independent  properties  in 
the  adiabatic  limit. 

The  growth  of  a  shear  band  and  the  dependence  of  its  width  on  material  properties,  such  as  strain 
hardening  and  thermal  conductivity  cbaracterbtics,  has  been  carried  out  by  Wu  and  FYeund  (1984). 
They  considered  an  imperfection-free  half-space  subjected  to  a  time  dependent  but  spatially  uniform 
tangential  loading.  The  half-space  material  exhibited  strain  hardening,  thermal  softening  and  strain  rate 


sensitivity.  They  applied  the  wave  trapping  idea  (cf.,  Erlich  et  al.  (1080))  in  the  numerical  simulation  of 
plastic  shear  wave  propagation  under  adiabatic  conditions,  as  well  as  with  heat  conduction  included.  It 
was  found  that  thermal  conduction  made  the  shear  band  slightly  narrower.  They  also  considered  both 
linear  dependence  and  lagarithmic  dependence  of  the  flow  stress  on  the  strain  rate,  and  found  that  the 
band  thickness  varied  with  rate  of  loading  differently  for  these  two  types  of  rate  sensitivity. 

The  mechanism  by  which  shear  bands  evolve  under  high  rates  of  loading  is  not  clear,  but  it  appears 
that  their  formation  is  quite  abrupt.  One  possibility  is  that  they  initiate  due  to  strain  concentrations  at 
material  or  geometrical  defects,  and  that  they  propagate  in  a  crack>like  fashion  to  become  the  bands 
observed  in  laboratory  specimens.  In  this  discussion,  the  phenomenon  is  studied  from  this  point  of  view. 
The  two  dimensional  antiplane  shear  deformation  of  a  large  block  of  material  is  considered.  As  a  first 
attempt  at  examining  the  shear  band  propagation  process,  the  material  is  assumed  to  be  hyperelastic 
and  incompressible  at  each  point.  Furthermore,  to  simulate  the  phenomenon  of  thermal  softening  due 
to  adiabatic  heating,  a  material  model  is  selected  which  shows  a  local  maximum  in  the  dependence  of 
shear  stress  on  shear  strain.  The  block  contains  a  small  region  in  which  the  stiffness  is  less  than  in  the 
remainder  of  the  otherwise  homogeneous  material,  and  this  small  region  is  called  the  defect.  Initially, 
the  block  is  supposed  to  have  no  strain  and  a  spatially  nnifwm  strain  rate.  At  points  remote  from  the 
defect  (compared  to  the  size  of  the  defect),  loads  are  applied  to  produce  a  constant  average  strain  rate 
for  later  times.  While  the  field  equations  do  not  appear  to  admit  a  complete  analytical  solution,  two 
simple  analyses  are  first  carried  out  to  shed  some  light  on  the  nature  of  the  transient  fields.  If  the  field 
equations  are  written  in  a  form  sutiable  for  analysis  of  plane  wave  propagation,  then  an  equation  for 
the  speeds  of  propagation  of  acceleration  waves  may  be  obtained.  The  direction  of  the  growing  shear 
band  is  determined  by  examining  the  inclination  of  an  acceleration  wave  which  becomes  stationary  as  a 
result  of  the  local  strain  field  attaining  the  critical  value  (corresponding  to  the  local  maximum  in  stress 
in  the  constitutive  relation).  In  addition,  the  field  equations  may  be  linearized  about  the  ongoing  uniform 
background  strain.  The  resulting  wave  equation,  which  is  linear  but  with  time  dependent  coefficients, 
contains  the  influence  of  the  defect  as  a  forcing  term.  The  wave  speed  resulting  from  this  wave  equation 
provides  an  estimate  of  the  speed  of  propagation  of  the  edge  of  the  shear  band  which  grows  from  the 
defect.  Finally,  a  complete  numerical  solution  of  the  problem  is  obtained  by  means  of  the  finite  element 


discretization  of  the  body  and  by  the  use  of  a  finite  difference  method  to  intergrate  the  differential 
equations  for  the  evolution  of  the  nodal  variables  in  the  finite  elements.  It  is  found  that  the  strain  in 
the  block  remains  more  or  less  uniform  up  untU  the  time  at  which  the  condition  for  instability  is  met  in 
the  defect.  Strain  is  rapidly  concentrated  at  the  edge  of  the  defect,  and  a  shear  band  grows  at  very  high 
speed  from  the  defect  out  into  the  remainder  of  the  block.  It  is  important  to  note  that  this  happens  for 
average  strain  levels  which  are  substantially  below  the  critical  strain  level  for  the  material  outside  of  the 
defect.  The  speed  of  propagation  of  the  edge  of  the  band  compares  very  well  with  the  estimate  of  speed 
obtained  from  the  linearized  analysis. 


vy 


1,  MATERIAL  MODEL 


In  the  linearized  theory  of  elasticity,  nontrivial  equilibrium  fields  of  antiplane  shear  are  possible  in 
the  absence  of  body  forces  for  any  homogeneous,  isotropic  material.  However,  according  to  the  discussion 
of  Adkins  (1954),  this  cannot  be  expected  for  the  exact  theory  of  antiplane  shear  deformation  in  finite 
elasticity.  In  contrast,  it  shows  that  there  must  be  some  restriction  on  the  strain  energy  function  of 
the  material  in  order  to  allow  for  nontrivial  two  dimensional  out-of>plane  displacement  fields.  Knowles 
(1076)  has  provided  the  necessary  and  sufficient  condition  for  a  homogeneous,  isotropic,  incompressible 
material  to  admit  nontrivial  states  of  antiplane  shear,  and  the  material  modeb  used  here  comply  with 
thu  conditon. 

The  second  restriction  on  the  material  considered  here  u  that  there  must  exbt  a  critical  strain  beyond 
which  the  field  equations  governing  equilibrium  fields  loose  ellipticity  for  homogeneous  deformations. 
Compliance  with  thu  restriction  will  allow  for  the  possibility  of  initiating  shear  bands  and,  for  antiplane 
shear,  the  existence  of  a  local  maximum  in  the  dependence  of  shear  stress  on  shear  strain  is  necessary 
and  sufficient. 

Following  Knowles  (1977),  we  consider  the  material  with  the  strain  energy  function  W{Ii)  defined 


where  A  u  the  first  invariant  of  the  left  Cauchy>Green  deformation  tensor.  If  u(zi,X2)  is  the  antiplane 
displacement  parallel  to  the  zs  direction,  then 


in  which  the  Greek  indices  take  the  values  1,2  and  repeated  indices  imply  summation  over  thu  range.  An 
index  following  a  comma  denotes  partial  differentiation  with  respect  to  the  corresponding  coordinate.  The 
material  description  (2.1)  involves  three  materia  paramaters,  namely,  the  shear  modulus  for  infinitesimal 
deformations  Ho,  the  hardening  parameter  n,  and  the  positive  dimensionless  parameter  6.  If  n  b  i  then 
the  material  b  linearly  elastic.  If  n  b  in  the  range  0.5  <  n  <  1  then  the  stress  mcreases  indefinitely 


with  increasing  strain.  And  if  0  <  n  <  0.5  then  the  stress  has  a  local  maximum  with  increasing  strain. 


For  the  ideal  material  defined  by  (2.1),  the  Cauchy  stress-strain  relation  for  antiplane  shear  defor¬ 
mation  is 

rs«  =  =  2W-(A)tt,« 

fas  —  2VV'(/i)u,„  u,a  -doiTi  +  tt)  -  di  (2.3) 

where  the  prime  denotes  differentiation  with  respect  to  A,  and  do  and  di  are  constants  to  be  determined 
by  the  boundary  conditions.  Introducing  the  quantities 


as  measures  of  the  shear  stress  magnitude  and  the  shear  strain  magnitude,  respectively,  in  the  body,  the 
constitutive  relation  between  r  and  k  for  the  material  being  considered  here  is 

+  **)  =  Hok  1-1- 

The  relationship  between  r  and  k  is  shown  for  several  values  of  n  and  (  in  Pig.  1.  As  can  be  seen  from 
the  figure,  r  is  a  monotically  increasing  function  of  k  tot  n  >  0.5.  On  the  other  hand,  there  exists  a 
local  maximum  in  stress  for  0  <  n  <  0.5.  The  value  of  1;  at  which  the  maximum  in  stress  occurs  is 


8.  FORMULATION 


A  rectangular  cartesian  coordinate  system  is  introduced  so  that  the  Xi ,  Z2-plane  b  the  plane  of 
deformation,  and  all  field  quantities  are  independent  of  xg.  In  addition,  there  exbts  a  region  around  the 
origin  of  coordinates  in  which  the  material  offers  less  resbtance  to  deformation  than  the  material  in  the 
remainder  of  the  body.  A  schematic  of  the  configuration  b  shown  in  Fig.  2.  Initially,  the  displacement 
b  zero  and  the  particle  velocity  varies  linearly  with  Xj  according  to 

w(0,x,,X2)=  axj  (3.1) 


Thus,  a  represents  the  uniform  average  strain  rate.  For  later  times,  the  velocity  field  at  points  remote 
from  the  defect  (compared  to  the  size  of  the  defect)  b  forced  to  have  the  form  (3.1).  For  points  in  the 
vicinity  of  the  defect,  however,  the  velocity  field  b  allowed  to  evolve  according  to  the  field  equations.  If 
the  body  were  homogeneous,  then  the  velocity  field  (3.1)  would  provide  a  solution  to  the  field  equations 
for  all  times.  Because  of  the  presence  of  the  defect  in  the  present  case,  strains  will  become  concentrated 
near  the  defect  and,  under  suitable  conditions,  a  shear  band  will  be  initiated. 

The  only  nontrivial  momentum  equation  for  antiplane  shear  b 


dri3  drti  _  5*u 


(3.2) 


in  which  t  denotes  time.  The  constant  do  in  (2.3)  b  equal  to  zero  because  the  end  faces  of  the  body, 
at  Xg  ^  ±L,  say,  are  free  of  normal  traction.  Substituting  from  (2.3)  for  the  stress  in  terms  of  the 
displacement  gradient  into  (3.2)  yields 


a 

dXa 


(3.3) 


which  b  a  second  order  quasilinear  partial  differential  equation  for  the  dbpalcement  as  a  function  of 
Xi,X2,f. 

One  of  several  ways  of  introducing  the  concept  of  characterbtic  surfaces  (cf.,  Courant  and  Hilbert, 
1962)  b  as  a  surface  ^’(zi ,  zg,  ()  =»  0  for  which  a  welFposed  Cauchy  problem  cannot  be  solved  in  general. 


When  this  is  the  case,  the  normal  derivative  d^u/d0^  will  be  indeterminate  on  this  surface.  Consequently, 
certain  fields  may  be  discontinuous  across  the  surface  without  violating  the  conditions  imposed  through 
the  field  equations. 
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If  new  independent  variables  are  introduced  so  that  the  V  coordinate  is  normal  to  a  prospective 


characteristic  curve  in  the  t  =  constant  plane,  and  the  other  "spatial”  coordinate  is  along  the  curve, 
then  the  equation  (3.3)  takes  the  form 


[2W*tl),a  i>,a  +4lV*'tt,a  +  f?  =  0 


The  p'lantity  R  represents  the  terms  which  can  be  determined  from  the  Cauchy  data.  FVom  its  definition, 
the  curve  V*  =  0  is  a  characteristic  curve  if  the  coefficient  of  in  (3.4)  is  equal  to  zero.  Thb 


condition  may  be  written  in  the  form 


(IV*  +  2W"u,?  )(^y  +  (IT  +  2W\l  )(^y 


+4W«u,i  tt,a 


di>  di> 


dxi  dzi 


If  the  direction  na  and  the  speed  V  are  now  introduced  as 


d^fdZg 


di>/dt 


then  Ra  is  obviously  the  unit  normal  to  the  spatial  characteristic  curve  0  »  0  and  V  b  the  corresponding 
speed  of  propagation.  The  curve  ^  =  0  m^  also  be  interpreted  as  an  acceleration  wave  (cf.,  Tfuesdell, 
1961).  Furthermore,  due  to  the  nature  of  antiplane  deformations,  the  speed  of  propagation  of  thb 
wave  does  not  depend  on  the  particle  velocity.  The  substitution  of  (3.6)  and  (3.7)  into  (3.5)  gives  the 


characterbtic  condition 


{W*  +  2VV«u,»  )nj  +  (»r  +  2W'u\  )n*  +  4iru„  u,*  runj  =»  pV* 


The  condition  for  the  formation  of  a  shear  band  could  be  viewed  as  being  equivalent  to  the  condition  for 
the  acceleration  wave  to  become  stationary,  that  is,  for  its  speed  to  become  zero.  Therefore,  it  is  possible 
to  find  the  direction  of  the  stationary  wave  front,  or  the  plane  on  which  a  band  forms,  by  setting  V  ssQ 
in  the  characteristic  condition  (3.8).  Hence, 


ni  2W"«,,  tt,2  + 

na  ”  2W"u,l 


(3.9) 


For  the  particular  form  of  the  strain  energy  function  being  used  here,  (3.9)  becomes 

2(1  -  n)^«,i  u,a  ±\/-{l  +  ^u,o,  u,*  )(1  +  (2n  -  l)^u,|,  u,/, ) 

—  = - - - - - IT - 1 - - -  (3-10) 

"*  (2n  -  l)-tt,f  +1  +  -tt,| 
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Inspection  of  (3.10)  and  comparison  with  the  stress>strain  relation  (2.5)  reveals  that  real  values  of  the 
ratio  fti/na  are  first  possible  for  values  of  strain  corresponding  to  the  maximum  in  the  r—k  relationship, 
that  is,  for  u,a  u,^  s=  n/(5(l  — 2n)].  This  implies  that  the  two  dimensional  problem  of  antiplane  shear  has 
basic  features  in  common  with  the  wave-trapping  concept  mentioned  in  the  Introduction  and  invoked 
by  Wu  and  FVeund  (1984).  Moreover,  it  can  be  demonstrated  that  the  equation  of  equilibrium  suffers  a 
loss  of  ellipticity  when  this  same  maximum  point  m  the  r  —  k  relationship  is  reached  (Zee  and  Sternberg, 
1983). 


In  addition,  it  is  evident  that  the  value  of  the  ratio  iti/nj  in  (3.10)  will  depend  on  the  boundary 
conditions  in  any  actual  deformation  field.  For  the  problem  at  hand,  the  background  deformation  field 
has  a  large  strain  component  in  the  xj  direction  and  no  strain  in  the  xi  direction.  Consequently,  it  can 
be  anticipated  that  ni/rts  =  0  and  the  shear  band  will  grow  in  the  X|  direction.  This  expectation  is 
borne  out  by  the  results  of  the  finite  element  calculation. 


One  special  case  of  interest  is  when  the  exponent  n  is  equal  to  1/2.  As  can  be  seen  from  (3.10),  the 
ratio  ni/fiQ  is  then 


^  ^  26u,i  tt,2  ±y/-l  -  2bu,a  u,tt 
ns  l'f26u,| 


which  is  always  imaginary  and,  therefore,  no  stationary  waves  can  exist.  Effectively,  the  same  conclusion 
can  be  reached  by  means  of  the  numerical  calculation,  and  this  is  further  discussed  in  section  5. 


4.  FINITE  ELEMENT  FORMULATION 


The  field  equations  are  satisfied  in  the  weak,  or  variational,  sense  of  finite  elements.  The  lengths 
At  and  A2  are  large  compared  to  the  size  of  the  defect.  The  boundary  conditions  are  chosen  so  that 
the  far-field  particle  velocity  (3.1)  is  imposed  on  Xi  =  h\  and  X2  ^  Aj,  the  symmetry  condition  of 
zero  displacement  is  enforced  on  rj  a  0  and  the  lymmetry  condition  of  zero  normal  derivative  of 
displacement  is  enforced  on  Z]  s  0. 

The  stress-strain  relation  of  the  material  and  the  equation  of  momentum  balance  are 

f.  * 

f4  =  /I  1  +  -tt,a  U,a 
L  ^ 

= pv,t 


subject  to  the  boundary  conditions 

v(t,zi,0)  —  0 
v{t,xi,hi)  —  ahi 

(4.2) 

V,l(f»0,Z2)  =  0 
v(t,Ai,Z2)  =  0 

The  initial  condition  is  given  by  (3.1).  In  the  foregoing  equations,  u  is  the  particle  displacement  and  v  is 
the  particle  velocity,  respectively,  in  the  Zs  direction.  Also,  in  these  equations  and  in  the  analyse  that 
follows.  Tat  is  replaced  by  fa  The  material  constmits  b,  n,  and  p  are  the  same  as  those  introduced  in  the 
preceding  section,  except  that  now  p  varies  with  position  zi,Z2  to  simulate  the  presence  of  the  material 
defect.  The  particular  form 


p  —  Po(l  +  «F(zi,Z2)i 


(4.3) 


is  assumed  for  this  variation,  where  po  is  the  constant  infinitesimal  shear  modulus  of  the  material  outside 
of  the  defect  and  the  function  F,  which  vanishes  outside  of  the  defect,  describes  the  material  variation 
within  the  defect.  The  shape  function  u  chosen  to  have  continuous  first  derivatives. 

It  is  convenient  to  introduce  nondimensional  coordinates  and  field  quantities  as  follows: 


’5' ^ *P  'yj '> 'Ji ’> ■?  ■> *> 'J  V •> ■> ’> ’JI ■■)*’>  V’> ■■'»,^’J -.T’T-g -^r^T^T^: 

m 


**< 


K<v 

\s 


E<v 


Eriy 


P:v? 


<1  ■»  —  (2  “  T-  >J  “  T-‘ 
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tv  cs  **  *  ® 

fcjTa  dtl  Colo 


(4.4) 


where  Cg  is  the  speed  of  inOnitesimal  elastic  waves, 


(4.5) 


and  p  is  material  mass  density.  For  n  <  0.5,  7g  is  the  critical  strain  at  which  the  shear  stress  (2.5) 
reaches  its  maximum  value.  In  terms  of  the  material  parameters,  70  is  given  by 


7* 


[i(l-2n)] 


1/2 


(4.6) 


If  stress  is  elimated  from  (4.1)  then  the  field  equations  may  be  written  in  the  form  of  a  first  order  system 
of  partial  differential  equations  in  terms  of  the  nondimensmnal  variables  as  follows. 


k-i 


(4.7) 


where  a  bar  under  a  subscript  is  used  to  indicate  that  it  is  to  be  viewed  as  an  exception  to  the  summation 
convention,  and  subscripts  following  a  comma  denote  differentiation  with  respect  to  the  corresponding 
nondimensional  coordinate.  The  boundary  conditions  are 


^il,6.0)-0 


Cola 


^(l|,l,^2)-0 


(4.8) 


in  terms  of  the  nondimensional  variables,  and  the  initial  conditions  are 


ah2 


Colo 


«'(0,(,,«8)-0 


(4.9) 


In  order  to  obtain  a  numerical  solution  of  the  equations  above,  Galerkin’s  principle  is  applied  to 
satisfy  the  equations  in  a  weak  sense.  Both  sides  of  the  first  equation  in  (4.7)  are  multiplied  by  a  function, 
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S>: 


say  which  is  arbitrary  except  that  it  must  have  the  same  smoothness  properties  and  must  satisfy 
the  same  kinematic  conditions  as  and  then  are  integrated  over  the  rectangular  region  in  Fig.  2.  This 
region  is  divided  into  a  40  by  40  array  of  like  rectangles,  each  of  which  is  composed  of  two  constant 
strain  triangles.  Application  of  the  divergence  theorem  leads  to 

j  dA  «=  -  y^(l  +  </)|l  +  ^(A27o)*«',o  A“®1  hlK^U),0  4^*0  dA  (4.10) 

where  /({i ,  (2)  is  the  imperfection  shape  functioa  viewed  as  a  function  of  the  nondimensional  spatial 
coordinates.  The  rectangular  region  is  then  divided  into  triangular  finite  elements  as  indicated  in  Fig. 
2,  and  both  the  displacement  field  and  velocity  field  are  expressed  approximately  in  terms  of  their  nodal 
values  in  the  usual  way,  that  is, 

(4.11) 

XO  B  V){0t 

where  the  interpolation  functions  6t  varies  linearly  within  each  element,  and  it  has  the  value  one  at  the 
t  —  th  node  and  the  value  lero  at  all  other  nodes.  Summation  over  all  nodes  in  the  finite  element  array 
is  assumed  in  (4.11).  Likewise, 

**  =  (4.12) 


In  view  of  the  fact  that  0*  is  arbitrary,  the  variational  statement  of  the  governing  equations  (4.10)  msy 
be  reduced  to  a  system  of  ordinary  differential  equations  in  t;  for  the  nodal  values  of  the  field  quantities, 


A/(rodr„  —  /?<(u>/:  e) 

®rie  ■■ 


(4.13) 


where  is  the  constant  (lumped)  mass  matrix.  The  system  of  differential  equations  (4.13)  is  solved 
numerically  by  means  of  the  forward  finite  difference  scheme  as  discussed  by  Zienkiewicz  (1977).  Numerical 
experiments  were  performed  to  establish  limits  on  the  time>step  range  required  for  stability  of  the  scheme 
and  to  test  the  accuracy  of  the  scheme  for  calculations  which  involve  many  time  steps. 

In  the  calculations,  a  form  for  the  imperfection  shape  function  is  required.  In  all  calculations,  it  was 
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,  ,  \iii<dandii<d 

/U».&)-{  (4.14) 

lO,  it  ii>  dor  d 

The  valve  d  »  0.1  was  used  in  generating  all  of  the  results  being  presented  in  this  report.  Thb  implies 
that  the  size  of  the  imperfection  is  approximately  one-tenth  of  the  size  of  the  region  over  which  the 
governing  equations  are  solved  numerically.  Because  the  shiq>e  function  (4.14)  has  a  continuous  dnivative, 
howevn,  the  material  weakness  is  actually  confined  to  a  somewhat  smaller  region.  Three  different  values 
for  the  strength  of  the  imperfection  t  were  selected,  namely,  c  «  —0.25,— 0.5,  and  —0.75,  for  calculations 
done  with  the  particular  value  of  hardening  (or  softening)  parameter  n  -■  0.4.  This  corresponds  to  the 
relative  rate  of  adiabatic  softening  observed  in  high  strength  steels  (cf.  Culver(lfi73)).It  is  found  that  the 
speed  of  propagation  of  the  band  is  not  sensitive  to  the  strength  of  the  defect,  but  the  strength  c  does 
infiuence  the  magnitude  of  the  strain  within  the  region  of  deformation  localization.  The  difference  is 
due  mainly  to  the  fact  that  a  shear  band  initiates  at  a  smaller  level  of  background  strain  for  a  stronger 
imperfection  (larger  magnitude  ct  c),  so  that  more  strain  has  accumulated  at  any  given  time.  The  results 
in  the  following  sections  are  baaed  almost  entirely  on  calculations  done  with  (  =  —0.5. 


S.  RESULTS 


In  order  to  get  a  better  nnderstanding  of  the  influence  of  the  parameter  n  on  the  propagation  speed 
of  a  shear  band,  calculations  were  performed  for  severrd  values  of  n,  namely,  n  »  0.4,  0.45  and  0.5. 
The  loading  rate  parameter  akifcario  wm  assumed  to  have  the  value  0.2.  Physical  parameters  which 
would  result  in  this  value  are,  for  example,  an  elastic  wave  speed  of  3000m  «~*,  an  imposed  strain  rate 
of  3000«~*,  a  transverse  dimension  of  the  block  0.2cm,  and  a  critical  strain  level  of  0.01. 

First,  the  distribution  of  particle  velocity  ^  and  shear  strain  'jtt  along  the  line  (i  «  2d  are  shown 
for  the  ease  n  a  0.4.  The  nodal  values  of  the  velocity  ^  are  shown  for  different  times  in  Fig.  3,  and 
the  variation  of  the  strain  component  728  (which  is  uniform  within  an  element)  is  shown  in  Fig.  4.  The 
stages  in  the  nucleation  of  a  band  are  evident  in  Fig.  3.  Initially,  the  particle  velocity  is  nearly  linear 
across  the  section.  As  the  deformation  advances  and  the  differences  in  material  stiffness  becomes  evident, 
material  particles  accelerate  for  points  near  the  defect  (due  to  reduced  resisting  stress).  The  diminishing 
resistance  of  the  material  near  the  defect  to  deformation  causes  an  unloading  wave  to  propagate  into 
the  block.  The  strain  dbtribution  shown  develops  accordingly,  as  shown  in  Fig.  4.  The  uniform  strain 
distribution  gives  way  to  a  distributioa  which  is  highly  localised  in  the  region  of  the  defect  and  which 
shows  relaxation  at  points  away  from  the  defect.  The  general  features  of  the  distribution  shown  in  Fig. 
4  are  similar  to  those  observed  in  the  one-dimensional  analysis  considered  by  Wn  and  Freund  (1984). 

In  order  to  investigate  the  propagation  of  the  shev  band  in  the  zi  direction,  the  distribution  of 
velocity  ^  in  the  (1  direction  for  (2  *■  i>  shown  in  Fig.  5  for  several  values  of  the  time  parameter  y. 
In  addition,  the  distribution  of  strain  723  in  the  (1  direction  for  the  first  row  of  elements  (0  <  (2  <  <1/^) 
is  shown  in  Fig.  6  for  several  values  of  the  time  parameter  y.  Fkom  Fig.  5,  it  is  evident  that  the  particle 
velocity  increases  rapidly  at  a  material  point  as  the  edge  of  the  band  approaches,  and  that  the  velocity 
decreases  equally  rapidly  after  the  edge  passes  by  the  observation  point.  The  variation  of  particle  velocity 
shown  in  Fig.  5  is  reminiscent  of  the  correspoinding  result  for  a  particle  near  the  fracture  plane  of  a 
growing  crack  as  the  crack  tip  passes  by.  The  r^id  variations  in  particle  velocity  suggest  that  material 
inertia  may  pity  a  significant  role  in  the  shear  band  growth  process. 

Because  they  were  unable  to  observe  partially  developed  shear  bands  in  their  experiments,  Costin 


et  al.  (1977)  concluded  that  the  speed  of  propagation  of  the  edge  of  the  bands  must  be  very  fast,  even 
on  the  time  scale  of  a  dynamic  experiment  conducted  with  the  Kokky>bar  apparatus.  To  obtain  an 
estimate  of  the  speed  of  propagation  of  the  edge  of  the  bands  for  the  present  model,  the  data  in  Fig. 
6  were  used.  It  was  noted  that,  after  the  band  was  wdl-developed,  fixed  levels  of  strain  propagated  at 
almost  constant  speed.  For  example,  if  attention  is  focused  on  the  strain  level  0.02,  it  is  found  that  the 
speed  of  propagation  is  approximately  0.57c;,  for  n  «  0.4.  If  the  same  calculation  is  done  for  n  »  0.45, 
then  the  speed  of  propagation  is  found  to  be  about  0.48Co.  Recall  that  Co  is  the  speed  of  infinitesimal 
elastic  waves. 

In  the  Appendix,  a  simple  analysis  is  described  in  which  the  governing  equations  for  the  problem 
at  hand  are  linearised  around  the  ongoing  back^onnd  strain.  It  is  found  that  the  displacement  field  b 
governed  by  a  linear  wave  equation  with  time  dependent  coefficients;  see  (A.3).  If  the  speed  of  propagation 
of  dbturbances  in  the  zi  dbection  b  estimated  from  the  equation  at  that  time  at  which  the  coefficient 
of  the  term  vanbhes,  then  the  foUowmg  result  b  obtained. 


(5.1) 


The  propagation  speed  obtained  from  thb  simple  approach  was  quite  close  to  the  speed  observed  from 
the  computed  fields  for  all  cases  tested,  always  within  10%.  For  example,  for  n  -■  0.4,  the  estimate  from 
(5.1)  b  0.584c,;  for  n  «  0.45,  the  estimate  b  0.517ew. 

In  order  to  get  a  view  of  the  complete  fields  at  a  fixed  time,  results  in  the  form  of  Figs.  7  and  8  were 
generated.  Fig.  7  b  intended  to  show  the  complete  partble  velocity  dbtribution  as  a  surface  over  the 
Cl,  Cs  plane  at  a  fixed  instant  of  time.  The  particle  velocity  at  a  point  is  the  elevation  of  the  surface  over 
that  point.  The  concentration  of  velocity  around  the  edge  of  the  band  b  evident  from  the  figure.  Fig.  8 
shows  the  dbtribution  of  the  strain  component  72s  for  the  same  time.  The  strain  is  clearly  quite  uniform 
over  the  whole  block,  except  for  points  on  the  shear  band  where  the  strain  levels  are  many  times  higher 
than  the  ongoing  background  strain. 

To  illustrate  the  complexity  of  the  stress  hbtory  experienced  by  a  material  particle  due  to  the 
passage  of  a  shear  band  edge,  the  triqectory  of  the  stress  point  in  stress  space  b  shown  in  Fig.  9  for  the 


eleineiit  adj^ent  to  the  symmetry  plane  {2  s  0  at  (i  a  2d.  In  the  early  stage  of  the  process,  the  stress 
component  ri  is  essentially  zero  and  the  component  rj  develops  within  the  nniform  ongoing  strain  field. 
As  the  band  edge  approaches,  the  stress  state  deviates  from  being  proportional  but,  some  time  after  the 
edge  of  the  band  has  passed,  the  state  again  becomes  one-dimensional.  Actually,  the  stress  hbtory  in 
the  second  row  of  elements  from  the  symmetry  line  (2=0  for  (1  «  2d  shows  a  greater  divergence,  as 
shown  in  Fig.  10.  This  stress  history  resembles  the  stress  hbtory  experienced  by  a  material  particle  near 
the  fracture  plane  due  to  passage  of  the  edge  of  a  crack. 

A  full  set  of  results  was  also  computed  for  n  0.45,  and  the  general  features  described  for  n  = 
0.4  above  were  again  obtained.  The  calculations  were  repeated  for  n  b  0.5  which  is  the  value  of  the 
hardening  parameter  for  which  the  stress*  strain  relation  has  a  horizontal  asymptote  for  very  large 
values  of  strain  but  the  stress  does  not  have  a  local  maximum  at  any  finite  strain  level.  In  thb  case, 
no  stationary  acceleration  wave  can  exbt  so  that,  based  on  the  approach  taken  in  section  3,  a  shear 
band  b  not  anticipated.  The  strain  distributions  analogous  to  Figs.  4  and  6  for  n  «  0.5  are  shown  in 
Figs.  11  and  12,  respectively.  Obviously,  the  material  defect  acts  as  a  strain  concentrator  and  the  strain 
intensifies  at  points  near  the  defect.  However,  a  shear  band  which  propagates  aw^  from  the  defect  b 
not  observed  in  thb  case.  Instead,  a  diffuse  nonlinear  zone  (like  a  crack  tip  plastic  zone  in  a  hardening 
material)  gradually  expands  away  from  the  defect. 


0.  DISCUSSION 


A  nonlinear  elastic  material  with  an  imperfection  under  finite  antiplane  derformation  has  been 
modeled  to  illustrate  the  process  of  formation  and  propagation  of  a  band  of  intense  shear  deformation. 
Several  observations  of  a  general  nature  can  be  made  concerning  this  analysis.  For  example,  the  effect 
of  the  imperfection  on  the  ongoing  background  deformation  field  is  very  small  until  the  strain  within 
the  defect  region  reaches  a  critical  value  for  instability.  The  resulting  load  transfer  to  material  adjacent 
to  the  defect  causes  a  strain  concentration  which  rapidly  forms  into  a  narrow  band  or  layer  which 
subsequently  propagates  in  a  crack-like  fashion  through  the  homogeneous  material  outside  the  defect 
region.  Because  the  material  has  no  viscosity,  the  thickness  of  the  layer  would  be  vanbhingly  small  if 
the  field  equations  were  satisfied  pointwise.  In  the  present  case,  with  the  resolution  of  the  analysis  set 
by  the  finite  element  mesh,  the  thickness  of  the  band  or  layer  is  always  found  to  be  one  mesh  spacing. 

Numerical  experiments  were  conducted  in  order  to  investigate  the  infiuence  of  the  strength  of  the 
initial  imperfection  on  the  evolution  of  the  shear  Ityer.  Of  conrse,  it  was  found  that  if  the  critical  strain 
for  instability  within  the  defect  were  increased  toward  the  critical  strain  for  instability  in  the  otherwise 
homogeneous  block  of  material,  then  the  shear  band  was  nucleated  at  a  higher  level  of  background  strain. 
However,  it  was  found  that  the  speed  of  propagation  of  the  edge  of  tbe  shear  band,  once  formed,  was  not 
sensitive  to  the  level  of  background  strain  at  which  the  band  formed.  That  is,  once  the  band  was  formed, 
its  growth  seemed  to  be  sustained  by  the  natural  strain  concentration  at  the  edge  of  the  band,  no  matter 
what  the  level  of  background  strain.  If  this  feature  carries  over  to  the  description  of  the  process  for 
more  realistic  material  modeb,  it  would  east  some  doubt  on  the  applicability  of  one  dimensional  shear 
models  in  the  study  of  band  formation  in  materiab.  The  one  dimensional  modeb  require  that  some  sort 
of  instability  criterion  be  met  globally  in  order  for  a  band  to  form,  whereas  the  two  dimensional  band 
growth  model  incorporates  a  natural  strain  concentrator  at  the  edge  of  the  band  so  that  an  instability 
criterion  must  be  met  only  locally  for  the  band  to  evolve. 

Finally,  the  limitation  of  the  present  anidysb  to  elastic  (that  b,  inviscid  and  path-independent) 
material  response  b  a  serious  deficiency  in  cmisidering  application  of  the  results  in  the  interpretation  of 
experimental  data  on  adiabatic  shear  band  formation  in  real  materiab.  It  b  anticipated  that  the  inclusion 


of  strain  rate  sensitivity  in  the  description  of  material  response  would  inhibit  the  formation  of  the  band 
at  the  defect,  whereas  the  inclusion  of  plasticity  effects  would  provide  stabilization  of  the  band  growth 
process,  analogous  to  its  influence  in  plastic  crack  growth  in  fracture  mechanics.  The  investigation  is 
continuing,  with  primary  attention  being  devoted  to  the  incorporation  of  more  realistic  material  models 
in  the  analysis. 
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APPENDIX 


If  the  strength  of  the  imperfection  t  is  small  compared  to  one,  then  a  regular  perturbation  analysu 
m^  be  carried  out  to  obtain  linear  differential  equations  for  the  deviation  of  the  field  quantities  from 
their  values  for  the  uniform  background  strain.  Due  to  the  presence  of  the  imperfection  described  in 
(3.4),  the  displacement  field  is  assumed  to  have  the  form 


u{t,Xi,X2)<^  (UZ2  +  tU{t,Xi,X2) 


{A.I) 


The  strain  components  corresponding  to  (A.1)  are 


U,i^€.U,l  u,2  =»  of +  ct/, 2 


(A.2) 


The  expressions  (A.I)  and  (A.2)  are  substituted  into  the  governing  equations  (3.2),  and  all  terms  are 
expanded  in  powers  of  e.  The  terms  of  order  sero  balance  identically.  A  necessary  condition  for  the 
terms  of  order  e  to  balance  is 


l/,n  +D U,22  —AU,tt  “  —aiF ,2 


(A.3) 


where  F  b  the  imperfection  sh^e  function  as  defined  in  (4.3).  The  coeficients  A  and  B  depend  on  the 
ongoing  background  strain,  which  implies  that  A  and  B  are  functions  of  time.  For  the  specific  material 
model  being  considered  here  they  have  the  form 
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(A.4) 


If  n  >  0.5  then  the  corresponding  homogeneous  equation  b  always  hyperbolic  because  B  is  always 
greater  than  aero  no  matter  how  large  the  uniform  background  strain.  Consequently,  all  plane  waves 
have  real  speeds  of  propagation.  However,  if  n  <  0.5  then  B  goes  to  lero  at  some  time  and  (A.3)  loses 
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hyperbolicity  at  this  time.  The  corresponding  wave  speeds  are  imaginary.  This  simple  argnment  suggests 
the  instability  of  deformation  which  results  from  attainmg  strain  levels  corresponding  to  the  part  of  the 
T—k  curve  with  negative  slope,  and  it  effectively  demonstrates  that  initially  small  disturbances  will  grow 
large  enough  to  significantly  perturb  the  homogeneous  background  strain  (see  abo  Hayes  and  Rivlin, 
1061)  and  give  rise  to  the  localized  shear  band. 


While  it  cannot  be  regarded  as  a  rigorous  step  in  apy  sense,  an  estimate  of  the  speed  of  propagation 
of  the  edge  of  the  band  may  be  obtained  by  examining  the  equation  (A.3)  at  the  instant  of  time  at  which 
B  »  0.  At  this  instant 


U,tt  ~^U,u  "  —'teF,i 


(A.5) 


where  A  >■  p/|/(0(l  +  The  "wave  speed”  for  propagation  in  the  Zi  direction  suggested  by 

this  equation  is  given  in  (5.1). 


FIGURE  CAPTIONS 


Figure  1:  Effective  shear  stress  r/|i«  vs.  effective  shear  strain  k. 


Figure  2:  The  plane  of  deformation  showing  a  coarse  representation  of  the  finite  element  mesh  and  the 
region  of  the  defect. 


Figure  3:  The  nondimensional  velocity  0.4  vs.  the  normalized  distance  (3  for  (1  »  2d  and 

several  values  of  normalized  tinm  9. 


Figure  4:  Shear  strain  'lat  for  n  »  0.4  vs.  normalized  distance  (3  for  (1  =  2i  and  several  values  of 
normalized  time  q. 


Figure  5:  Nondimensional  velocity  ^  torn  ■■  0.4  vs.  (1  for  (3  b  d/4. 


Figure  6:  Shear  strain  ^at  for  n  «  0.4  vs.  distance  (1  at  (3  «  d/8. 


Figure  7:  The  velocity  surface  di  the  {t>{z  plane  at  time  17 «  3.6.  The  Imear  velocity  distribution 
in  (3.1)  has  been  subtracted  out. 


Figure  8:  The  shear  strain  surface  73$  over  the  (1,  plane  at  time  7  ■■  3.6. 


Figure  9:  Time  history  of  stress  in  the  stess  plane  at  (1  ■■  2d  and  0  <  (3  <  d/4  for  n  »  0.4. 


Figure  10:  Time  history  of  stress  in  the  stress  plane  at  (i  2d  and  d/4  <  (3  <  d/2  for  n  «  0.4. 


Figure  11:  Shear  strain  733  for  n  *■  0.5  vs.  distance  (a  »t  (1  ■■  2d  for  several  values  of  normalised  time 

n- 


Figure  12:  Shear  strain  733  for  n  »  0.5  vs.  distance  (1  at  (3  ■>  d/8  for  several  values  of  normalised  time 
V- 
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